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Many physical structures are threatened by combinations 
of loads of varying magnitudes from various sources. For 
instance, bridges, piers, dams, and buildings can experience 
loads frem wind, snow, ice, tides, earthgquakes, and so 
forth. In many instances the total load or stress exferi- 
encec Fy a structure varies in time in an apparently randon 
Pasi ten Certain components of the loads vary rather 
Slowly; others occur more nearly aS impulses, such as those 
associated with winds or earthquakes. The problem is to 
design a structure to withstand the superposition of random 
loads- from many sotrces with at least an approximately 
understocd probability. 

The purpose of this thesis is to describe and investi- 
gate certain Simple but somewhat realistic probabilistic 
loai models for use in design, and perhars safety assessment 
of structures. In this thesis we confine attention tc the 
Superposition of just two load types; shock loads, and 
constant loads. For example, winds gusts, flash floois, and 
earthquakes have varying magnitudes and have relatively 
short durations in ccmparison to the times between their 
eccurrences. These will be modelled as instantaneous shock 
loads. Cn the other hand, snow, ice, or water accumulation, 
or even the presence of slowly moving vehicles present loads 
that remain nearly censtant “in time, occasiona iy semana 
to new levels. These will be modeled as constant loads that 
change infrequently. Throughout this investigations weaeE 
Le assured that the effective stress exerted by several 
types of loads acting simultaneously can be expressed as a 
linear comkination cf the component loads, the load 


components being treated as stochastic processes. 
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The times between changes in magnitude of the constant 
1oad frocess are independently and identically distrikuted 
with distribution function H. The succeSSive nagnitudes of 
the constant load are independently and identically distrib- 
Meed with distributicn function F. Given the constant load 
PIOCESS, the shock load process iS a compound Poisson 
process. The successive shock ioad magnitudes are indepen- 
dently and adentically distributed with distribution func- 
won. GC; The conditicnal rate of arrival of shocks, given 
that tte magnitude cf the constant load process iS x, is 

AMX) - 

Let X(t) be the magnitude of constant load process at 
time t; and Y(t) be the magnitude of shock load process at 
time t; Z(t)=X(t)+Y(t) is the superposition of the two loads 
at time t, and i {t)=s0F Z(S), the maximum load combination 
Piece t: |, See Ligure ted. 

Let T,={inf t20:; Z(t)>x } ie, the first-passage time 
Sor the loac combination to exceed a stress level x. Ty 
will represent the time to faiiure of a structure whose 
strength is x, and is subjected to a stress history 
c2 (t eee 2 0} . 

Gaver and Jacobs (1981) studied the above nodel for the 
case in which the distribution H is exponential and the 
shock load rate A’ 1S aconstant independent of the shock 
10ad process. Related load combination models that have 
been studied include Peir and Cornell (1973) (Ref. 1] Wen 
(1977) { Ref. 2] Pearce and Wen (1983) [Ref. 3]. 

Asyrrtotic results which appear in Gaver and Jacobs 
(1984). and Jacobs (1984) will be described in chapter 2 for 
the distribution of Ix aS x3 and for the tail of the 
ete Ibutron Of tf, for finite x. 

im “Chapter 3, a simulation program for the load 


combination model will be described. 


mm mr mm mm me ee ee ee ee 








ERS de) SSS a i OS Si ee ER coc! og a ee ES ee I ee i I egy Sm Ln a 
— 
a. 
See 
SS cer) A es ee) ete) a, A OS pees EY geet EY ree) CE OD ee Gig Red ol pees cme cata & a? Yee 
. 


Pagqurem lid The Load Combination Process. 
In chapter 4, a Simulation study of the use of the 


asymptotic results Ofe chapter 2 to approximate the 


distribution of Ty, will be described. 
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A. THE DISTRIBUTION CF Ty, FOR LARGE X 


This section summarizes some of the results found in 
Gaver anc Jacobs (1981) [Ref. 4] for the model in which the 
Gistribution H is exfonential(Aa) and the shock load arrival 
rate A is constant independent of the constant load process. 

sue eye 4 we consider the distribution function of the 


(aii Ooaascomolnam@ on co occur during (0,t1, M(t). 


H(t) =P (M(t) <x} 
=P{M(t) <x, Tr >t} +P (M(t) <x,T <t}o ree (2-1) 
where T; iS the arrival time of the first shock. 
A renewal argument yields; 
P (M(t) <x,T, >tp= E> exp (-ue Gex-y)) Ecdy) 
Pim (tj) <x,T, <3 =f rxedu|* exP(-av Glx-4)) Hy t-v) ECAw) 
where G(x)=1-G (x). 


Next take Laplace transforms with respect to t of (2-1). 


Hele) = fo enFt Hace) de 


SITE og CE, NSE ON GE a alla a (222), 
where 
eee Sealine iG) x-men (GY) =a =~ ~~~ = (2-3) 
neues 
Sy = _ Me (8B) a te ET Qi 
memseetiomea be daftreult te invent H,(¢) for any anter- 
ee wenomce Or the Gdistributiors F(x) and G(x). However, 


useZul information can still be gleaned from equation (2-4). 


Tinswme bom Cie aetTipfatvion of iI, , 


EQM(t)<xJ=P(Ty>t) 000 tee o eee eee een e----------- (2-5) 
chus, Bo@Ws POS (Gh Sie) Gis (2-6) 


Demme s 20, ih edilation (2-6) to f1rnd that 
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< My (0) 
= hore = mm mmrwrwmnm rrr nn wr wee = 
m (x) B (Ty ) H,(0) iy -~A Mx(0) (2 up) 
The next result is the limiting property for the first 
time the load combination process exceeds a given level x. 
The limiting distribution of Ty is exponential in the sense 


faa 


xk 


a) scale 4m =e. OO” eee eee - 
a F{m(x) T, >t} =e (2=a3 
here x, =1nfi fe, *G(t)— a 


the preef of this result for a more general model that 


. 


* 


includes that of chapter 2 can be found in Jacobs (1984) 
[Reme 5 |e 


Be THE TAIL OF THE DiStTRIBUTION Client EN) Baw 


In this section we will summarize results for the asymp- 
totic behavior of P(I,>t) for large t and fixed x as t7?@ 
for the model in which the distribution H is exponential and 
the shock@foad arrival rate 1S Gonseant: More details and 
results for a more oaeneral model that includes that of 
Chapter 2 can be foune@™in Dacobs (19384). 

From the equations (2-1) and (2-5) 

P(T, >t) = e-rtf* Ecdyy e- HEHE 


trfr ers ds [*ecdy) eMGU-HS pee ytes) 9 ----- (BS) 
Let 
L(ds)={* Fcdyy ner ars ds == (2-10) 
r . Bee 
Le {? Fd) eazy io @ ae g)1t} —------- oe 1h 
r» 
bis fe6) | 8 Ecdy) Ue ene mee 1) 
The renewal equation (2-9) gives; 
P (Be >t) g(t) + ico syria (2-1) 
where g(t) = pcdy) @TAtHGr$Jt = = (2- 14) 


Following the deveiopment in Feller (1971; page 376) 
fRef. 611t will be assumed that there exists Janu ere. 


such that 
4 ee exw4 1 ¢d4) Sn a el ee er (215) 


flee 


Deeseeroct isSweimaque and Since the distribution of L is 


defective, iS positive. 


Let 
ft vagqQae’: . fey > eee aac 
Coe acon eee ---- 3175 
S22. See eae 
Then 
p* (T, >t) =g*(t)+ [PE (7, > (t-s)}L¥ (ds) 9 --------- (2-19) 


It 1S assumed that the assumption for the key renewal 


theorem holds. Afrlication of the key renewal theorem 


yields 
: Koot Pea ee ee oe 2 
Lin e P(T, >t) = alg (t)at (2-20) 


(g#(t) at= f F (dy) Meee a TAt Ma (x-$) ) t at 


mex (ESTE ee = 
=f F(dy) ATK Glx-4) - K(x) eee) 


Rewriting the equation (2-15), the equation determining 
K Ais 2 

) x xX = - 
1= { ES | Etdyy > eT LATHE U-DIS Ye 


 — aa 
~ Fo A + AGLX-Y) — KY 
and 


We ace aXe eA Ne AB) dt 





= ie ECAY) 22 ae 
2 [A+ A a@tx-$) - Ko }* 


the key renewal theorem implies; 


/ 


Lim e Kat D le Pee = a a oe ee ee oe oe ats (Tis 24) 
a Ni) Se a 
Dd = 
mean ples. [A+ UM G(X-4) - KOW)* 


It appears to he difficult to solve for Kk analytically 
in general. We will discuss the numerical computation. of K 
for two examples. In both examples, N=A=1 and the 
distribution orf the shock and constant load magnitudes F and 
Gare of the forn; 

Eitey = een 


(Oe ears 


ee 


In example 1, a=r=1 and in example 2, a/b=2. 


TO Compute the kK for F(x) = feuaecee— eX P 
let 
A 
ex es eee. ee eee ee as 
E(x) =i Fag) KN + AL GOX-$) - KY WZ) 


where A=A =1. 
Fer model 1, 
f (x)= TRE UP 0 7)-2 "Ln 1-40) = Ain (1-0) 2% e% 3 -- (2-26) 
for K<l+e" jane voy. 
ror K=1, £(X) haS a Singular point, so by uSing the 
ft HOPITAL aun Ret ee 
Zin £(x)= Lin g'{K) 7h (x) 
where g@) is the numerator and h(x) is the denominator 
ef the sunctivonm =i 
For kK =1, 
f (1) == (e -é*) ---- (2-27) 
on aes” ZG 
iQue ss as [Lucxeerby? _ fuser e)* | ie Sayre i) 





+ OT On (pee 07hX) — Lm (Uezne + e-bx)}} ----- (2-28) 


for «<1te™ and K #1. 
For K=1 £() has also singular point, so by using the same 


rmethcd as used in model 1; 


for K =1 

A ere a 
Differentiation of the ecuation (2-25) gives 

=a(K) =) eae eee (2-30) 


[ A+ a eP&B -KJ> 
EGUatiON (2-30) 41s > g¢eatemmetnanme tos O<x<1+e% and K*¥1. So 
the function f(x) is a “Monotone increasing function in the 
interval 0<x<1l+e™ “) WSing this) west we | compumomtnce% 
from the equation (2-22) by applying the bisection method 
[Ref. 8}. 
Using the computed k we compute the constant y¥* such that; 


oo nN a P| ee a er 
he ae C A+ GK-$) -KJZ 


1S the ccnstant of equation sizo24 
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Fer case a=b=1, 


y¥s ea ine haere”) - Sar] 
—-\ o-x =V y= -x oe a oD 5 ae ee mecca 
fe-* (2-K)-2e-* In (E™(2-K)) UGE }} ( ) 


For case a/b=2, 


4 2 n —bx 2 =e ~bx ~bXx ai 
i er a ~ fucker ene —3e*(uu-xycr-e-**)) 


3 Z / A 
$3072 *Cim (1 eee) Ln (UK ee OO ge aee rave y\ 27 32) 





ITI. COMPUTER SIMULATION 


This section describes the computer simulation model. 


we denote the constant load arrival rate bya, the shock 


load arrival rate by gr, “the distribution of comearanc toad 


magnitude by F(x), and the distribution of shock load magni- 


tude by G(x). The simulation model will be used to study 


the quality of the arproxmamations of the disthrburionyore 


by the asymptotic results described in chapter 2. 


The following describes the simulation algorithn. The 


stress ievels are X, <xX,<X,...-.---- ee 


Tl, 
2, 


Set 2=0, No =U. 

Set V=0, generate a constant load magnitude,C, and 

duration interval Tc. Find the largest index i1>N, 

Siche tha wear 

PUt Nec iia re. tic tee cee enc cl 

set 
ly =f [esor pie ie, 

Se Myth) - 

Generate a shock load magnitude,S, and a time until 

the shock occurs,Ts. 

If the shock interval does not exceed the excess life 

of the constant load duration,Tc, find the largest 

index 2>N, Suc beth awee Crs. 

Let N, be that index. 

Set . Ty: =TEVtT Ss Weecr <r ,. 

set No =N,\ ,V=Vtis, [C=fests., Go etoustep 3. 

If the shock interval exceeds the excess life time 

of the constant load duration, . 


set T=T+Te¢, go to steum.. 
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One realization cf the Simulation is completed when one 


Ty Poecon ea CamHol scocnem Level x, , |\Si<n. ‘Hack sirulation 
consisted of 5000 replications. Further details are 





G 


ae eae, ee) 


| 


gee Seen | 


ee 


Figure 3.1 Generating T. 


iveiudeds 1m “chapter 5. 
The key subprograms are following. Other programs will 


be explained by reading from the output results. 


RANDC ; generate the Ix. 

CCMPUT ; compute the P(Tx=0), true P(Tx=0), sample mean, 
true mean, standard deviation, coefficient of 
variance. 

Sample prebability 
F (Tx =9) = ese is obtained from 
the sample. 

Sr uempreotabalit y 
E (Tx =0) =1-P(Tx>0) =1-P (M(0) $x) =1- JEP (ay) =F (t) 1 
obtained from the equation (2-5). 

The sample mean 
a 


mn 
=— Ee 
X ne, TX; and 


a7 


the true mean E(T,) is obtained in special 
cases from the equation (2-7). 
See details in Gaver and Jacobs (1981). 

Sample standard deviation 


“a 


nN 
Ss i TAeae) is obtained from the sample, 





and 

Sample coefficient of variance 
COD) a is obtained from the sample 
mean and Sample standard deviation. 

CCMP ; compute tke sample guantiles, and the guantiles 
cf an exponential distribution with mean E(Tyx). 

A sample quantile is 
9 = (dxsanple Size) th order statistic of sample and 
fa = -E(te) dm ( $=) is the guantile of exponential 
distribution week mean 5.0) 

CCM ; compute the quantile of the sample conditional 
distribution given Tx>0O and the exponential distrib- 
ution with mean E(T,)/P(Ty>0). 

The sample guantile of the conditional distribut- 
10n 
9 e=[ Ad x (number O£ positive Ty) |th erdenm st arrs rie 
is obtained for sample cf poSitive T,. 

The conditional quantile of the exponential dist- 
Tibution with mean E(Tx)/P(Tx>0) is 


a BCiOss al 
bac = Tea ee ke) = 


SK 3 conpute the Kesvleneenat. 
K={ K(x) 2 JP Fedy) 
bisection method. 


A * 
Nesbit eene 
INT ; compute the constant y* such that ; 


ae Me cd ee eee 
: JoF a CA+ M&M Gx-$) -K jz 
by substituting thesk et oundarnes.. 
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the following cases were simulated. 








aUSiteLDution Sasc cS AACC) F (xX) G (x) 
a,exvonential ie 1 1 1-e” 1-6" 
(a=b=1) Tae 1 . iT 1! 
; ? 1 LA La 
b,éexponential ie 1 1 1-e 7* ete Gane 
(a=2, b=1) Ly 1 2 . i 
Se ? 1 1" 11 
Gpex cenentaa | 1 1 wie 1-e~* 1-e™% 
We brea yang 25 1 en sf i 
arrival rate 3, 1 eS Ny i 
d,Fareto case tS 1 2 1- =x ies 
2. 1 2 eae 1-e™% 
‘ 
oP 1 2 i Sarre 


These different settings will be referred to as case (A-1), 
oe nee o= lt) 7 ember) yee (B-3), (C-1), (C-2), (C-3), 
eee = 7 slo oy respectively throughout this paper. for 
each above case, we ccmpare the asymptotic result (2-8) and 
the experimental results. 

The next setting is to examine the exponential approxi- 
Pe OmeonnEnemm tai. Cr the daStribution of T, for finite xX. 
The cases of identical exponential distributions, 
Poole i—ce vanasdifterent exoonential distributions, 


F (x)=e™ G(x)=e” where a/b=2 


x 


were considered. 
The following conditions are considered. 
é; A =A =1 x) Gin) —en* 


ae DS ST Piel cac Gx) ee 


se 


These different settings will be referred to as case (£) and 
case (F) respectively. For each case, compute the #&, and 
constant %§* and compare the guantiles of the asymrftotic 


results (2-24) to the Simulated data for each level of x. 
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IV. ANALYSIS OF THE SIMULATION RESOLTS 

Bo Hditcatredamam chapter 2 the exacta@astribution of T, 
Pwcmioetclu lt tomontain in general. In this chapter sinula- 
tion will he used to study the accuracy of the approxima- 
moms Of the G@stribution @er Ty by the exponential 
distributions suggested by the asymptotic results of chapter 


<.@ 


A. THE DISTRIBUTION CF Tx FOR INCREASING X 


the limiting result (2-8) indicated that the distribu- 
tion of Tx approaches the exponential as level x increases. 
This result suggests that the distribution of Ty can be 
approximated by an exponential distribution at least for 
large x. In this section this exponential approximation will 


te studied via a simulation. 


In the first collection of three models to be 
SOEs lderca ms (x) —G (x) - €*. The times between constant load 
changes are exponential with parameterAand the shock ioad 
arrival rate mw iS a constant independent of the constant 
idad process. 
For these cases it is possible to derive an analyt- 
ical expression for E(Tx). 
ee. Nee) ~Uxe® +e Im (RE) | 
ae [4+ Hx - = Im SS ? 
In the Tables, P(Ty, =0) refers to the simulated 


See memonamiity |that Tx =0> X-BAR Gives. the simulated 


sample mean of Ty, ; ST-DEV is the simulated standard 
deviation of T,; VAR(X-B) tS. the variance cametic Sanere 
mean; COEF-V is the sample standard deviation divided by the 
Sample wean. T P(Txy=0) 1S the theoreticai probability that 
T.=0, F(X); #T x-bar is the theoretical nean of Tx computed 
from the equation (4-1). 

The Simulated coefficient of variation in all three 
cases decreases as the level x increases becoming guite 
close to the theoretical exponential value of 1 when x=5. 

To compare the simulated distribution of Tx to the 
approximating exponential distribution, quantiles from the 
Simulated data were computed and compared to the corre- 
Sponcing quantiles cf an exponential distribution having 
theoretical mean (4-1). The Simulated quantiles were 
computed as described in Chawiter 3a The exponential 
-d-qguantile is given Ey 

gq, =~a-ln (1-4) 
where a is the mean of the exvonential. 

For each level of x, the first row in the tables 
(Awdjiy (A-2), (AnS)) Lommthe euantales of theses enon tachi 
ix gives the quantiles of the Simulated data and the second 
row gives the corresponding exponential quantiles for an 
exponential distribution having theoretical mean (4-1). 

Cne way the distribution of Tx differs from an exfo- 


nential is that it has an atom at Ol inaeparciculan: 


Di(ty 0) =) oe 
for the case considered here. 

Quantiles were used to compare the simulated condi- 
tional distribution of Ty given Ty820)) to ane “pemen aa 
daseralutionaiavaing mean (Tx )/P(Tx>0). The simulated ‘opel 
tiles for the conditional ~~ drstri1ourron mewerecrconedc-1a 
Cescribed in chapter 3. These cuantiles were compared to 
those of an exponential distribution having theoretical mean 


GE yas (ye pale) 


Ze 


i eeicwt diabes (AS) (Ae ziy andeGosy) for the guan- 
Pos Or tiemecOnGimtemeial. Ji Stiabutionsitor each level x, the 
first row shows the simulated quantile and the second row 
the corresponding approximating exponential guantile. 

The guantiles of the Simulated conditional distribu- 
tion are much closer to their exponential approximation than 
the guantiles for the unconditional distribution to their 
exponential approximation. However the quantiles for the 
unconditional distritution get closer to those of their 
approximating exponential as x increases. 

Gonparangeeea bles WA), (Ne2), and 9(A-3), it appears 
that increasing the arrival rate of shock or the rate of 
Change of the constant load magnitude decreases the gquan- 
tiles. The change of arrival rate of shocks appears to have 


mac greater effect. 


2. cOnStant and sSkock Load Maqnitudes Have Different 


In the next three cases, tie distryputicon we: tie 
shock load magnitudes is exponential with parameter 1 and 
the Gistrikbution of the constant load magnitudes 1s exponen- 
tial with mean 90.5. The other model assumptions in tables 
(Bs Vy (a2), (B.3) correspond to those in tables (A.1), 
faveez2 yp ame (As3) . 

As before, for each level x, theweias tAacONean the 
quantile table gives the simulated guantile and the second 
row gives the approximating exponential quantile using the 
exponential distribution with theoretical mean. Comparing 
the guantiles of in table (B.1) with (A. 1) (respectively 
(emeeereeen (fees 6Ccande(ee))06C With ) «6©9©(f%.3)), it 1s seen that 
Gecreasing the mean constant load magnitudes has increased 
PiomGucnemlesm—NErimti rer, 2f dppears that the convergence of 
MiemeGtetribtwimonme or ix tO am exponential is faster in the 


case of smaller mean constant load magnitude. 


De) 


The exponential approximation to the conditional 
distribution of Ty given Tx>0 has theoretical mean 
E(Ty ) ABeGIy > 0). Once again for small levels x,it appears 
that the exponential approximation to the = ¢omdi tionan 
distribution of %is better than that for the uncendi tiomam 


One. 


—— ee eee ee ee = ss os oe — —_— [= SSS = = — —_— a =—_ — = —_— 


load Maqnitude 


In the next 3 cases, constant load magnitudes and 
Shock load magnitudes are exponential with mean 1. Constant 
ioad magnitudes change according to a Poisson process with 
rate 1. Given the constant load magnitude at time t is C, 
the prokakbility a shock load will arrive in the time 
interval [t2tth ) sec) bos oni. 


In case (C-1) the ccnditional shock load arrival 


Fate ws ue tee Comparison with table (A.1) indicates that 
conditional shock arrival rate wic)" =C has the exftectmor 
cecréeaSing the quantiles of fTy,. Further the exponential 


approximation to the quantiles of Tx, 1S not as good as in 
case (A-1). 
The approximating exponential distribution to the 
conditional distribution of Tx given T,>0O has a mean of 
E (te) Feetly 20). The quantiles of the exponential approxina- 
tion to the guantiles of thebcornditional disteibutieneaprcan 
to be closer than those for the unconditional distrikution. 
In case (C-Z) snock loads arrive with conditional 
acrival rate Ade)" =expiee Comparing the table (Geajmesaed 
(C.Z), it 1s seen that the quantiles of the Case Gime are 
less than those for the case u(Cy=exp(C). This is because 
interarrivail times tends to be larger for the case 
H(G\ =e one) Comparing tables (A.1) and (C.2) indicates 
that the exponential approximation to the quantiles of 7x is 


not as good as in the case M=1. 
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ie case wea) thescondittenal shock arrival rate is 

W(Cy' =exp (-C). Comparing tables (C.2) and (C.3) indicates 
that the guantiles fcr the case 

A(C)t =exfZ (-C) are less than those for the case y(C)' =exp (C) 

and furthermore the exponential approximation of Tx appears 

Bomee ClLOSC@REOE the case (C-3) than that of for the case 

fe-2). In all of the cases the simulated mean was used as 


the parameter in the exponential approximations. 


In the next three cases, the times between constant 
load changes of magnitude are exponential with mean 1. Shock 
loads arrive according to a Poisson process with rate 2. In 
case (D-1) the constant load magnitudes have a Pareto 
distribution with parameter 4=1 and the shock load nmagni- 
tudes are exponentiaily distributed with mean 1. Comparison 
with table (A.2) indicates that the Pareto constant loa4 
Cee idence ye=tect of decreasing the gwantiles of Tx. 
Further the exponential approximation to the gquantiles of t& 
1S not good in the Pareto case; (the approximating exponen- 
tial has a mean of tke simulated sample mean). 

The approximating exponential distribution of & 
given T, >0 to the conditional distribution has a mean of 
Poy )7P (Ty 29). Once again the guantiies of exponential 
approximation to the guantiles of the conditional distriku- 
tion appear to be closer than those for the unconditional 


G@usthie eat lore 


im case (D-2) the constant load magnitudes havea 
Pareto distribution with parameter A=2. Comparing takies 
(eat ce 2 eetteis seen that the quantiles of the simu- 


lated distribution cf Ty for the case d=2 are larger than 
trose for the case c=1. Comparing tables (D.2) and (A.2) 
indicate that the exponential approximation to the quantiles 
Cou ives mumeote as good for the Pareto case as for the 


exponential case. 


ZS 


In case (D=32) the constant load magnitudes have 
Pareto distribution with parameter “«=1 and the shock load 
magnitudes have Pareto distribution with ¢=1. In all three 
cases the simulated mean was used in the approximating gquan- 


tiles. 


Bo THE EXPONENTIAL TATE OF THE DPESTRIBUTION OF Tren 
FINITE X 


In this section simulaticn will be used to study the 
exponential approximation 


suggested by the asymptotic result 
tim -PCTer’) = 


{30 Ties aie 
This 1S an approximation for P (1, >t) Mer fixed finite eee 
St:iould become more accurate as t becomes large. 

Two caseS are considered. In both cases shock loads 
arrive accordjiing toa Poisson process with rate 1 and 
constant load magnitudes change at the times of arrival of a 
Poisson process with rate 1. The shock load magnitudes have 
exponential distribution with mean 1. In case (EP, the 
constant loai magnitude distribution 1s exponential with 
mean 1; In case (F) it 1s exponential with mean 0.5. 

Lescription of the computation of *« and T¥* for these two 
cases can be found in chapters 2,3. It can be seen from the 
computed value of x(x) anda F(x) appearing in tables, (Exnt9 
and {F.1) that K(x) 2S avppoaeching thesvalue (7: Gg (where 


E(Ty) 1s the theoreticai mean) as x becomes larger and for 
/ 


—— eee y(x) 1s approaching 
E(%) 


the cases computed x(x)< 
1 as x becomes larger. 

-To study the exponential approximation to the distribu- 
tion of TIT , guantiles from the simulated data were computed. 
These can be found in tables (£.2) and (F.2). For each level 
x, the first row presents the simulated quantile; the second 
row gives the approximating é-quantile computed by 


Tg yeh bm Uae) os) 


A 


The third row shows the approximating &-quantile computed by 
Ga =-E(™)In(1-&) 0 === == --------------------- (T ) 
where the theoretical mean is used for E(T,). 

Comparing the guantiles in tables (E.2) and (F.2) it can 
ke seen that the appfroximating guantiles (LZ) are close to 
the Simulated ones even for 2.40. The approximating quan- 
tiles computed in (zg) do less well but improve as x becomes 
larger. The two approximating guantiles approach one ancther 
as x becomes larger, as expected. Approximating guantile (rz) 
has the advantage,however, of being eaSier to compute. 

Table (G.1) presents 90 % confidence intervals for 
A-~guantiles with A420.9 for case (E). These confidence inter- 
vals were computed using a large Sample approximation, see 
GComever (T9307 page 111-114) [ Ref. 9}. 

Comparison with table (E.2) suggests that the approximating 
guantiles approximates the Simulated ones well for these 


values. 


Cia GONCLUSTONS 


Approximating guantile ZL iS aiways better than approxi- 
Tada iuatieruesiu ALPEOxXiMation — approximates q, well for 
aS Smail as 0.40. However ~~ 1S more difficult to compute 
thaw e. 

The performance of approximating quantile ZI can be 
Ih PEOVECeDyeeGmancing it toe 


gy = _ ECT) Son (1 - ) 


PC Tx >) 
and using it to approximate the quantiles of the conditional 
CistmeraoutiOoneoum sy given t0>0; for the models considered 


PimESOmeiseecasyetorvcempute being equal to F(x). 
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V. SIMULATION ALGORITHM 


VARIABLE DEFINITICN 


Det ee Os a(t) ek} 

where Z(t)=superposition of constant load and shock loai 
Peace UD etee mele ve TMX: e NECK CL <. oye wee ees <* 
Ncol=maximum number of repetition of generating Ty 
Lc=constant load magnitude 

LS=sheck load magnitude 

Tc=constant load interval 

Ts=shock load interval 

Tpce=ccunting the number of repetition 


L=counting the maximum level of X; 


AIGOFRITHM 


Infut;parameters (zr,m, a, b, seed numbers ) 
ieee aero mo rmeoo revel X; 5 X CX.cX,<....-.205- OGRE: 
i1pc=0 
Fepeat 
set T= 
M=0 
TTs=0 
Repeat 
generate ccnstant load;Lc 
generate constant interval;Tc 
Pawideeesicn that L={Sup 1;Lc2x; } 
Set Tx, = T0 for 1S1isL 
M=L 
Ci=Tc 
A, generate the shock load;lLs 
generate the shock interval;Ts 
Eee -TS>Ca, then 


TO] 04 Ie 

TSs=TSs-C1 

TTs=0 

generate constant load;Lc 
generate constant interval;Tc 


Find 1] such that h={sup 2: tec 
> 


Tf LOM, then 
Set Ty: —=2) VPOreeic re 
M=L 
Else continue 
Ci=Te 
GO 1LO.E 
Else, Fand lesuen tind & pe SDE SOLOS hay 


TTs=1Tis+Ts 
Set Tx, =T.+TTs for M<aisL 
Ci=Cieus 
Go to A 
Jntil (L=Nrow) 
Ipe=iper 
Until (Ipc=Ncol) 


End Aigorithn 
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